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Abstract. A strong negative dependence property for measures on {0, 1}" - 
stability - was recently developed in [5], by considering the zero set of the 
probability generating function. We extend this property to the more general 
setting of reaction-diffusion processes and collections of independent Markov 
chains. In one dimension the generalized stability property is now indepen- 
dently interesting, and we characterize the birth-death chains preserving it. 
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1. Introduction 

In statistical physics a fundamental object of concern is the partition function, 
with its zeros having special relevance. For example, by introducing the effect of 
an external field, the partition function becomes a polynomial in the external field 
variable. As exemplified by the Lee- Yang circle theorem in the case of the Ising 
model [12], the general location of partition function zeros can indicate possible 
phase transitions. 

A related object in probability is the probability generating function. However, 
the locations of its zeros were little studied before the recent work of Borcea, 
Branden, and Liggett. In [5j, a strong negative dependence theory for measures on 
{0, 1}™ was obtained; in particular, it was shown that if the generating function 
(in n variables) has no zeros with all imaginary parts positive, then the measure is 
negatively correlated in a variety of senses: negative association, ultra-log-concave 
rank sequence, Rayleigh property, and others. 

The classification of linear transformations preserving the set of multivariate 
polynomials that are non-vanishing in circular regions was recently resolved in 
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[2], with the investigation providing a general account of such polynomials and 
unifying several Lee- Yang- type theorems [3]. 

Our results are as follows: Using this framework we will generalize the neg- 
ative dependence result in [5] to measures on {0,1,2,...} S - S countable - with 
application to independent Markov chains and reaction-diffusion processes. The 
one-coordinate case is also independently interesting; more specifically, the proba- 
bility measures under consideration can be decomposed into a sum of independent 
Bernoulli and Poisson random variables. 

Call such measures on {0, 1,2,...} t-stable. (The formal definition is given 
in Section 3.) In the last section we characterize most birth-and-death chains 
preserving this class of measures: 

Theorem 1.1. The birth-death chain {X t ;t > 0} with an infinite number of non- 
zero rates (e.g. irreducible on N) preserves the class of t-stable measures if and 
only if the birth rates are constant and the death rates satisfy 5 k = d\k + c^fc 2 for 
some constants d\ 1 d2- 

One example is the pure death chain with rates 5k = k(k — l)/2, which 
expresses the number of ancestral genealogies in Kingman's coalescent - a well- 
studied model in mathematical biology [TTJ [101 El]- In particular, by taking the 
initial number of particles to infinity, we obtain that the number of ancestors at 
any fixed time has the distribution of a sum of independent Bernoulli and Poisson 
random variables. 

2. Stability and Negative Association 

We first review the relationship - established in [5] - between negative association 
and the zero set of generating functions for measures on {0, 1}™. 

Definition 2.1. A polynomial /(x) € C[x] = C[xi, . . . ,x n ] is called stable if / ^ 
on the set 

H" = {(an, ...,x n ) G C" : Im{xj) > Vj}. 

Let © [x] be the set of all stable polynomials in the variable x. 

If / has only real coefficients, it is also called real stable. The corresponding 
set of real stable polynomials is denoted ©r [x] . 

Note that a univariate real stable polynomial can only have real zeros. 
One key fact from complex analysis is the (multivariate) Hurwitz's theorem 
on zeros of analytic functions: (see footnote 3 in [7]) 

Theorem 2.2. Let Q be a connected open subset ofC n . Suppose the analytic func- 
tions {fk} converge uniformly on compact subsets offl (normal convergence in the 
vocabulary of complex analysis). If each fk has no zeros in Q then their limit f is 
either identically zero, or has no zeros in O. In particular, a normal limit of stable 
polynomials with bounded degree is either stable or 0. 
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For /i a probability measure on {0, 1}™, let 
l 

f li (x 1 ,...,x n )= J2 tih,...,i n )x?—xZ =Wxf l) --.x$ n K (2.1) 

The last expression is just compact notation for the middle sum - the T)(i) are 
the coordinate variables for fi. f ^ is known as the probability generating function 
for fi. With this identification between measures and polynomials, we will freely 
abuse notation by referring to measures with stable generating functions as stable 
measures (such measures are also termed Strongly Rayleigh [5 j , by their connection 
with the Rayleigh property). 

The concept of stability easily generalizes to countably many coordinates - 
a measure \i on {0, 1} S is stable if every projection of fx onto finite subsets of 
coordinates is stable. 

While the definition of stability is purely analytic, it implies two strong prob- 
abilistic conditions. Recall that a probability measure fi is negatively associated 
(NA) if, for all increasing continuous functions F, G depending on disjoint sets of 
coordinates, 

J FGdfi < J Fdfi J Gdp. 

The following was proved in [5]: 

Theorem 2.3. Suppose / M is stable. Then /i is NA. 

The second (and less difficult) probabilistic consequence of stability is the 
following [T61I25]. 

Theorem 2.4. Suppose fi is a measure on {0, 1} S such that is stable. Then for 
any T C S, 

£■>(<) = 

ieT iGT 

where the Q are independent Bernoulli variables and the equality is in distribution. 

For S finite this latter result has been known since the work of Levy [14] - 
See also Pitman [18] for more combinatorial and probabilistic properties of stable 
generating functions, and the connection with Polya frequency sequences. 

3. Stable measures on {0, 1, 2, . . . } s 

Suppose n is a measure on {0, 1,2,... }™. The generating function of \x is now the 
formal power series 

oo 

fn(x 1) ...,x n )= ^2 (j,(ii,...,i n )xf ■■■x%. (3.1) 

ii,...,i n = 

If /i has finite support, then / M is a polynomial. In this case, let N be the maximum 
degree of / M in any of the variables x\, . . . ,x n . We will want to represent / M by 



4 



Thomas M. Liggett and Alexander Vandenberg-Rodes 



a multi-afBne polynomial. To do this, we recall the k-th elementary symmetric 
polynomial in m variables 

eo = l, &k{x\i . . . , x m ) := ^ Xi 1 Xi 2 • • ■ Xi k . (3-2) 

l<ii <i2 < - <ik<m 

Then for a univariate polynomial f{x) = XfcLo a kX k we define its N-th polarization 
as 

E /N\^ 

Pol N f(x x ,. . .,x N ) := 2J ^ k J a k e k (x!,. . .,x N ). 

The iV-th polarization of a multivariate polynomial is then defined to be the 
composition of polarizations in each variable. By considering x\ = • • • = Xn = x, 
notice that if PoImI € <8[x] then f E &[x]. The converse also holds: 

Theorem 3.1. (Grace-Walsh-Szego) . Suppose f has degree at most N. Then f is 
stable iff Pol is stable. 

Many proofs of this result and its equivalent forms exist; see the appendix of 
[3J, or [2U chapter 5]. 

Definition 3.2. We say that a function /(x) defined on C™ is transcendental stable, 
or t-stable, if there exist stable polynomials {/ m (x)} such that / m — >• / uniformly 
on all compact subsets of C™ (/ can then be expressed as an absolutely convergent 
power series on C"). Let (5[x] be the set of all t-stable functions - this is also 
known as the Laguerre-Polya class [13]. Let ®r[x] be the set of all real t-stable 
functions. 

We will again abuse notation and say that a measure \i on N n is transcendental 
stable (or t-stable) if its generating function lies in ©r [x] . Similarly, if /x has finite 
support and its generating polynomial is stable then we say that fi is stable. Of 
course, a stable measure is automatically t-stable. 

The papers by Borcea and Branden [21 14] characterized the linear transfor- 
mations preserving stable polynomials by establishing a bijection between linear 
transformations preserving n-variable stability and t-stable powers series in 2n 
variables. We will not require their full result here; however, the following char- 
acterization of t-stable powers series - the technical cornerstone upon which the 
above bijection rests - will be most useful. 

Recall the standard partial order on N n : a < (3 if on < Pi for all 1 < i < n. 
Then for a, j3 € N" and letting /3! = ft! • • • ft!, define 

(/^ = if^ \\ iia <^ = otherwise. (3.3) 

IP - «)' 

Theorem 3.3 (Theorem 6.1 of 2 ). Let /(x) = X^aeN™ c a xQ ^ e a formal power 
series in x with coefficients in K. Set /3 m = (m, m, . . . , m) G W l . Then /(x) £ 
C5r[x] if and only if 

/ m (x) := V (/3 m ) Q c Q e C5 R [x] U {0}, (3.4) 

a</3 m 
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for all m £ N. In this case, the polynomials / m (x) — > /(x) uniformly on compact 
sets. 

This classification also has the following immediate consequence: 

Corollary 3.4. The class 25r[x] is closed under convergence of coefficients. In par- 
ticular, the set of t-stable probability measures on N n is closed under weak conver- 
gence. 

Proof. Suppose that for each n, 

with Ca^ — > c a for each a. Then for each m the stable polynomials fm\~x), defined 
in (|3.4[) above, converge normally to the polynomial f m likewise obtained from /. 
Hurwitz's Theorem implies that each f m is stable, and applying Theorem l3.3l again 
we conclude that / £ 25 R [x] . □ 

By the following proposition, we can say that a measure on {0, 1, 2, . . . } s is 
t-stable if every projection onto a finite subset of coordinates is a t-stable measure. 

Proposition 3.5. The class of t-stable measures is closed under projections onto 
subsets of coordinates. 

Proof. It suffices to consider projections of n coordinates onto n — 1 coordinates. 
Suppose that /i is a t-stable measure on N". If f(x\, . . . , x n ) is its generating 
function, notice that the generating function of the projection of ijl onto N n_1 is 
f(x±, . . . , £ n _i, 1). By Theorem 13.31 it suffices to show that the approximating 
polynomials f m {x\, . . . ,x n -\,l) are stable. But this follows by considering the 
(complex) stable polynomials f m (xi, . . . , x n —i, 1 + i/k) and applying Hurwitz's 
theorem as k — > oo. □ 

We can now give an extension of Theorem 12.31 

Theorem 3.6. Suppose [i is a t-stable probability measure on {0, 1,2,... } . Then 
fi is NA. 

Proof. By the previous proposition and a limiting argument it is sufficient to 
show the result for measures on N ra . Let /(x) be the generating function of /i. By 
definition, / £ ©r[x]. Let {/a?(x)} be the stable polynomials converging to / as in 
Theorem l3.3l which we can assume are normalized so that /at(1) = 1- Let /j,n be the 
respective probability measures on {0, 1,2,..., N} n . Hence by the Grace- Walsh- 
Szego (GWS) theorem, Pol^fN is the generating function for a stable measure 
/2at on {0, l} nN . Let 

{Qf,l<i<ri,0<j <N} 

be the coordinates of fiN, such that 7^ = J^j dj is the i-th coordinate of ^m- Hence 
for bounded increasing functions F and G on {0, 1,2,...}™, depending on disjoint 
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sets of coordinates, we have 
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The inequality above follows because F(in + • • ■ + ieijvi 
an increasing function in the nN variables (similarly with G), and the Qj are all 
negatively associated by Theorem [231 The normal convergence of Jn — ► / implies 
the weak convergence /ijy — > fJ>, concluding the proof. □ 

We can also characterize all t-stable measures on one coordinate. 

Proposition 3.7. A probability measure on {0, 1,2,...} is transcendental stable if 
and only if it has the same distribution as a (possibly infinite) sum of independent 
Bernoulli random variables and a Poisson random variable. 

Proof. Suppose / is a t-stable generating function for a non-negative, integer val- 
ued random variable. By Theorem l3.31 / is a normal limit of univariate polynomials 
with all zeros on the negative real axis. An appeal to the classical theory of entire 
functions (e.g. [121 VIII, Theorem 1]) indicates that / can be expressed as the 
following infinite product: 



f(x) = Cx«e™ Y[ 



k=l 



ak 



for some q G N, a > 0, < 0, and ^2 \ a k\ 1 < oo. A little rearrangement - using 
the fact that f(l) = 1 - gives the following alternative expression: 



= x q e a{x - 



^ll[0--Pk)+xp k ] 



k=l 



where p k = 1/(1 — ctfc)- This we recognize as the generating function for the sum 
of a non-negative constant q, independent Poisson(cr) and Bernoulli (p^) random 
variables. Conversely, any generating function of this form with ^2pk < oo is 
automatically t-stable, as e x is the normal limit of the polynomials (1 + x/n) n . □ 

By projecting onto finite subsets of coordinates (taking limits if need be) and 
setting all variables in the resulting generating function to be equal, we obtain the 
following extension of Theorem 12.41 

Corollary 3.8. Suppose fi is a t-stable measure on N s . Then for any T C S the 

number of particles located in T - according to /j, - has the distribution of a sum 
of independent Bernoulli and Poisson random variables. 
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3.1. Markov processes and stability 

Suppose {rjt;t > 0} is a Markov process on N". We define the associated linear 
operator T t on power series with bounded coefficients by letting T t (x Q ) be the 
generating function of {rytl^o = a} for each a € N", and extending by linearity. 
This is well-defined because J2 k>0 P(rjt = k) = 1. 

Definition 3.9 (Preservation of stability). We say that a Markov process {i] t ; t > 0} 
on N" preserves stability if for any stable initial distribution, the distribution at 
any later time is t-stable. That is, the associated linear operator T t maps the set 
of stable polynomials with non-negative coefficients into the set of t-stable power 
series. 

The process rj t preserves t-stability if for any t-stable initial distribution, the 
distribution at a later time is again t-stable. That is, T t maps the set of t-stable 
power series with non-negative coefficients into itself. 

In fact, these two definitions are equivalent. 

Proposition 3.10. A Markov process preserves t-stability if and only if it preserves 
stability. 

Proof. Only one direction needs proof. Assume the process preserves stability. Let 
T t be the associated linear operator, and / = J2 a c a xa be the generating function 
of a t-stable distribution; hence c a > for all a and ^ Q c a = 1. By Theorem 13.31 
there are stable polynomials /„ = J^a c « xa with ci™^ — > c Q , all > 0, and 
with J2 a c « < 1- Suppose that 



Since probability is conserved, d Qjl g = 1, and hence by dominated convergence, 
^ c4 n) (i Q ,/3 — > ^2 c^a.fi as n -> oo. 

a a 

In other words, the coefficients of T t f n (which is t-stable by assumption), converge 
to the coefficients of T t f . T t f is then t-stable by Corollarv l3.4l □ 

We now give a couple examples. 

3.2. Independent Markov chains 

Suppose {X t (l), X t (2), . . . } is a collection of independent Markov chains on S with 
identical jump rates. Set 

i>i 

so that the resulting process is a collection of particles on S jumping independently 
with the same rates. This is well defined as long as rj t (x) < oo for all igS, t > 
- one possibility is to restrict initial configurations to the space E defined below 
for reaction-diffusion processes. 
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Proposition 3.11. The process {f]t',t > 0} preserves t-stability. Hence, assuming 
that the initial distribution is t-stable, the distribution at any time is negatively 
associated by Theorem \3.b\ 



Proof . Let fit be the distribution of rjt, with /io t-stable. We need to show that for 
each finite T C S, the projection /i t \T is t-stable. Taking finite T C Si C 62 C ■ ■ • 
with each S n finite and S n /" S, we can approximate \x t \r by the sequence /Uj |t, 

(n) 

with each /zj the distribution of the independent Markov chain process on S n 
given initial distribution ^o\s n and jumps restricted to staying inside S n . Hence 
by Corollary 13.41 we can assume finite S. 

Suppose now that S = [n], and only jumps from site 1 to site 2 are allowed. 
In this case, assuming a jump rate (7(1,2), each particle at x independently has 
probability p := 1 — e - * 9 ' 1 ' 2 ' of moving to y. Hence the associated linear operator 
T t takes 

x^x* 2 ■ ■ ■ x%> 1— > (px 2 + (1 - p)x 1 ) ai x^ 2 

that is, 

T t f(xi, . . . ,x n ) = fipx2 + (1 -p)xi,x 2 , ■ ■ ■ ,x n ), 

which preserves the class of stable polynomials. By permuting variables, this ar- 
gument also holds for general i,j e [n]. 

Recalling the Banach space Co(N n ) of functions that vanish at infinity, con- 
sider the strongly continuous contraction semi-groups S 1 '^^) on Co(N") defined 

by 

5«(i)/( 77 )=EV(^ i ), 

where ;t > 0} is the (Feller) process of independent Markov chains which only 
allow jumps from site i to site j. Then for a t-stable intial distribution /x, we just 
showed that /j,S l,J (t) - the distribution of rf^ assuming initial distribution /j, - is 
again t-stable. By Trotter's product theorem [9j p. 33], the process allowing the 
jumps {it— > j},{k > 1} has semigroup 

S(t)= lim \s^(t)s k > l (-)] n . 

Hence by Corollary 13. 4[ fiS(t) is again t-stable. Including all the possible jumps 
one-by-one, we conclude that the whole process {rjt ; t > 0} on N n preserves t- 
stability. □ 

Remark 3.12. Let A4 be the set of probability measures on N described by random 
configurations rj with coordinates 



n(k) = J2 My,- 



where li, Y2, ■ • • are independent random variables with values in S U {00}. In 
[T5] it was shown that A4 is preserved by the process of independent Markov 
chains, and that measures in M. are NA. Proposition 13.111 is a generalization of 
this result, since it is easily checked that the class Ai is contained in the class of 
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t-stable measures. Indeed, if /i G Ai and S = [n], then /i has a generating function 
of the form 

= l[[P( Y i = oo) + P(Yi = l)xi +-.. + P(Yi = n)x n ], 

i 

by the independence of the Y^s. Furthermore, (non-constant) product measures 
for which each coordinate is a sum of independent Bernoulli and Poisson measures 
arc t-stable, but are not contained in the class M.. 

3.3. Reaction-diffusion processes 

In addition to having the motion of particles following independent Markov chains, 
we can also allow particles to undergo a reaction at each site. Let p(i,j) be transi- 
tion probabilities for a Markov chain on S. Given a state rj € N s , we consider the 
following evolution: 

1. at rate /3V* a particle is created at site i, 

2. at rate 5 l r -, a particle at site i dies. 

3. at rate r/(i)p(i,j) a particle at site i jumps to site j. 

The most common example is the polynomial model of order m, where the birth- 
death rates for each site are 

m— 1 m 

p k =Yl W - 1) • • • (* - i + 1), s k = J2 d J k ( k - 1) ...(*- i + 1). 

3=0 3=1 

Reaction-diffusion processes originated as a model for chemical reactions [2D] , and 
subsequent work by probabilists has focused on the ergodic properties [HI El H] ■ 

To construct the process, we require a strictly positive sequence fcj on the 
index set S, and a positive constant M such that 

^~]p(i, j)kj < Mki, i € S. 
j 

Furthermore, the birth rates must satisfy 

^/^fc* < oo. 

i 

Then we take the state space of the process to be 

E = {rj eN s :^2v(i)h < oo}. 

i 

See [5J chapter 13.2] for the details of the construction. 

In a very simple case we have preservation of t-stability: 

Proposition 3.13. Suppose the reaction- diffusion process on N s is well-constructed 
with f3 l k = b l and 8 l k — d l k - this is the polynomial model of order 1 and site-varying 
reaction rates. Then the process preserves t-stability, and hence - assuming a t- 
stable initial configuration - its distribution at any time is negatively associated. 
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Proof. The strategy here is the same as with Proposition 13.111 To reduce to a 
reaction-diffusion process on a finite number of sites, we approximate using the 
construction in [fH Theorem 13.8]. Furthermore, on the locally compact space N™, 
the reaction-diffusion process with at most constant birth and linear death rates is 
now a Feller process, as can be seen from [9j Theorem 3.1, Ch. 8]. Hence by Trot- 
ter's product formula and Corollary 13.41 we only need to show that the following 
processes preserve stability on N": 

1. constant birth rate V at a single site i. 

2. linear death rates at a single site i (5 l k — d l k), 

3. jumps from site % to site j at rate r](i)p(i,j) 

For (1), we note that with a constant birth rate, at time t a Poisson(6 l i) 
number of particles has been added to the system - i.e. the original generating 
function is multiplied by e 6 **^- 1 ), preserving t-stability. 

(2) can be thought of as the process in which each particle at site i dies 
independently at rate d 1 . Hence the associated linear transform is defined by 

' ' ' i ' ' ' "^n ) — ' ' ' [1 6 ~\~ 6 '' X n 

As the affine transformation x t— > ax + (1 — a), (a > 0) maps the upper half plane 
onto itself, T t preserves preserves the class of stable polynomials. 

Finally, we note that (3) was already seen to preserve stability from the proof 
of Proposition [3~TT1 for independent Markov chains. □ 

Other reaction-diffusion processes do not preserve stability in general. On one 
coordinate, a reaction-diffusion process is just a birth-death chain, so by Theorem 
11.11 the only possible generalization would be to quadratic death rates. In this 
case, unfortunately, the associated linear transformation T t will not preserve stable 
polynomials with positive roots; indeed, assuming death rates 8k = k{k — 1), 
quadratic polynomials with a double root inside the interval (0, 1) will not be 
mapped to stable polynomials under T t . A much more complicated example - 
which we shall not reproduce here - shows that quadratic death rates on multiple 
sites does not preserve the class of stable probability measures. 



4. Birth-Death Chains 

Our goal in this section is to prove Theorem ll.il As just noted above, in the case 
of quadratic death rates the associated linear transformation does not preserve 
all polynomials with real zeros, and hence we cannot rely on the classification 
theory in [2\. From the probabilistic point of view, it would be useful to have a 
similar theory for linear transformations on polynomials with positive coefficients; 
however, in what follows we will make do with several perturbation arguments - 
the main idea being that a polynomial's roots move continuously under changes 
to its coefficients. 
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Recall that a (continuous-time) birth-death chain {X t ;t > 0} is a Markov 
process on the non-negative integers with transitions 

k i — y k ~\~ 1 at rate /3fc, k i— > fc — 1 at rate 5^, with £o = 0. 

We only consider rates {13k, 5k} such that the process does not blow up in 
finite time - see Chapter 2 of [T7] for the necessary and sufficient conditions, as 
well as for the construction of the process. 

Because of our definition of "preserving stability" (Definition I3.9|) . Theorem 
11.11 and the auxiliary results below only make sense for birth-death chains with 
infinitely many non-zero rates (e.g. when the chain is irreducible on N). Indeed, if 
one fixes n and considers rates ftk = n—k and 5k = k for < k < n, with all other 
rates zero, then it can be seen that this process preserves all stable polynomials of 
degree at most n. 

The generating function at time t is given by 

oo 
k=0 

By Theorem 2.14 of [T7], the transition probabilities 
p t (j,k) = P(X t = k\X =j) 
are continuously differentiable in t and satisfy the Kolmogorov backward equations 

j t Pt{j, k) = f3 jPt {j + l,k) + 5 jPt {j - 1, k) - {p s + 6j)ptV, k). 
By Theorem 2.13 of [T7], \p t (j,k) -p s (j,k)\ < 1 -p\ t - a \U,j)- It follows that 

r. OO , 

fe=0 

for \z\ < 1, provided that ^(0) is bounded. Iterating this argument, one sees that 
4>{t,x) is C 2 on [0, oo) x {z : \z\ < 1} if X(Q) is bounded. 

Proposition 4.1. Suppose the birth-death process {X t ;t > 0} preserves stability. 
Then there exist constants bo, b\, &2, d\, c?2 so that 

fik = bo + b\k + &2^ 2 and 5k = d\k + c^fc 2 - 

Proof. Suppose 

n 

<j>(Q,z) = c Y[{z- Zk), 
fc=i 

where c > 0, n > 3, and —1 < z\,...,z n < 0. By assumption, the generating 
function <fi(t, z) of X t has only real roots for t > 0. If z\ — Z2 — w is a root of 
4>(0,z) of multiplicity exactly two, and e is small enough that \z k — w\ > e for 
k > 3, then Rouche's Theorem implies that for sufficiently small t > 0, <f>(t, z) has 
exactly two roots in the disk {z : \z — w\ < e}. Therefore, for small t > 0, there 
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exist real z(t) so that (f>(t,z(t)) = and lim^o z(i) — w. By Taylor's Theorem, 
there exist s(t) £ [0,t] and y(t) between z(t) and w so that 

=4>(t,z(t)) = t^(0,w) + l -t 2 ^(s{t),y{t)) 

+ t{z{t) - w) ^L {s{tiv{t)) + \ {z{t) w ) 2 S (s(t) ' m) - 

Dividing by t and letting t 4- leads to 

d(j), , ld 2 (b, . , (z(f)-») 2 
2-^(0, w) (0,io) = -lim^-^ - < 0. 

Noting that 

H(0,u;)=2cn(*-a), 

fc=3 

we see that 

^(0,H (4.1) 

changes sign when Z3 crosses u>, and hence is zero when Z3 = w. 

To exploit this fact, we need to compute (|4.1[) . Recall the fc" 1 elementary 
symmetric polynomials e k {x\, . . . ,x n ) defined in (|3.2j) . If /1 is the distribution of 
Xq and A"o < n, then 

n n n 

0(0, 0) = £ = c]J(z- z k ) = C ^(-l) fc e fc (zi, z n )z n -\ 

k=0 fe=l fc=0 

so 

fi(k) = c(-l) n ~ k e n - k (zi, z„). 
Therefore for \z\ < 1, by the Backward equations, 

^(0, z) = f>(l - + + - M(l)(/3/ + Si)}* 1 

1=0 

n 

=c(l - z) ^(-l)"-'e n . t ( Zll z„)[4^" 1 - 

It follows that the expression on the right is zero if z = Z\ = Z2 = z% = w for any 
values of w, Z4, z n € (—1,0). In this case, 

e k (zi,...,z n ) = w l ( 3 J efc-j(z4, z n ), 

i ^ ' 

where i ranges from max(0, fc + 3 — n) to min(£:, 3), so 

^e n _ fc _ i (z 4 ,..,z n )[4^ fc+i - 1 -/W^] =0. 
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Interchanging the order of summation and letting k >->• k — i, we see that the 
coefficient of each of the e„„fc's is zero: 



(-!)%_<-&_<«>] = 0, 



or equivalently 4 - 34+1 + 34+2 - 4+3 = and 4 - 34+1 + 34+2 - 4+3 = 0. 
So the birth rates 4 an d death rates 4 are quadratic functions of k, and the 
boundary term 4 — enforces 4 = d\k + dik 2 for some di,d,2- □ 

With the next proposition we resolve the "only if" part of Theorem 11.11 

Proposition 4.2. The birth-death chain preserves stability only if the birth rate is 
constant. 

Proof. Assuming that the chain preserves stability, we will show the birth rates 
/4 satisfy 4 > 4+i f° r each fc, so by Proposition ^. II 4 is constant. 

By iterating the Kolmogorov backward equations, one can obtain the follow- 
ing approximations for small t > 0: 

t 2 

p t (k, k + l) = i(4 + o(l)), Pt (k, k + 2) = -(44+i + o(l)). 

Similarly, 

p t {k,k-j) = 0(l)t j , 

where O(l) denotes a uniformly bounded quantity, and o(l) — > 0, as t — > 0. 

Suppose that we start the chain with k particles; the initial distribution has 
generating function f(x) = x k . We also can assume that 4)4+1 > 0- Then the 
generating function for small t > will be: 

f t (x) = • • ■ + (1 + o(l))x k + (4 + o(l))tx k+1 + (44+i + o(l))^x k + 2 + ■■■ 

Since ft(x) is t-stable, by Theorem 13.31 the following polynomial has all real, neg- 
ative roots: 

- ■ + a + ■>(!))*' + + afegew . 

As the hidden coefficients are o(l), Rouche's Theorem implies that k roots of ft.k+2 
are also o(l). Thus the remaining two roots a, b satisfy 

, , 2(fc + 2)+o(l) 



ab — 



4+i* 

(fc + 2) 2 + o(l) 



44+it 2 

Solving for real a, 6 implies that the discriminant 

it~ 2 (k + 2) 2 [/?-^ - (44+i) _1 + o(l)] > 0, for small t > 0. 
Taking t — > 0, we conclude that 4 > 4+i- D 
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We now concentrate on the "if" part of Theorem 11.11 Here is a useful fact 
about quadratic death rates: 

Proposition 4.3. The generating function of the birth-death chain with rates /3k = . 
8k = k(k — 1) satisfies the Wright- Fisher partial differential equation: 

d d 2 

— 4>(t, z) = z(l - z)—(j)(t, z). 

Proof. Let p t (l,k) — P(X t = k\Xo = I) be the transition probabilities. Since the 
process stays bounded we can use the Kolmogorov Forward equation : 

j t Pt(l, k) = (k + l)k Pt (l, k+l)-k(k- l) Pt (l, k) (4.2) 

Then, since Ht{k) = J2i M(0P*(^ we nave 

d d 
—<f>{t,z)=^n{l)- Pt {l > k)z k 

k,l 

= KOK* + l)*Pt(i, k + 1) - k(k - l)p t (l, k)]z k 
k,l 

=z(l ~z)J2 Kl)Hk - l)pt(i, k)z k - 2 = z(l - z)-^4>(t, z). 
k,l 

□ 

Proposition 4.4. The birth-death chain with rates Pk — 0, 5k — k(k — 1) preserves 
stability. 

Proof. As above, let <f>(t, z) be the generating function of the chain at time t, and 
now assume that </>(0, z) has only real roots. Setting 

r = inf{t > 0; <j>(t, z) is not stable}, 

by Hurwitz's Theorem (Theorem I2.2p <p(r, z) is stable. Hence it suffices to prove 
that for any stable initial distribution there exists an e > such that <f){t, z) is 
stable for all < t < e. 

For any root w < of 4>(0, i), we can write 

0(0, «) = (*- «;)"<?(*)> ( 4 -3) 

where n is the multiplicity and q(w) ^ 0. We will show that for all small enough t 
the generating function (p(t, z) has n distinct real roots of distance approximately 
\ft from w. This follows directly from showing that 

<f>(t, w + at 1 ' 2 ) = t n ' 2 p{a) + o(t n / 2 ), as U 0, (4.4) 

where p(ct) has only simple real roots. Indeed, (|4.4[) immediately implies that for 
all small enough t, <p(t, z) changes sign n times near w. Since we can do this for 
each w, this shows that cf>(t, z) has only real roots. 
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Set k — [n/2\ . To show (|4.4I) . we Taylor expand 4>(t, z) in the first variable, 
i.e., for small t > 



^(m) = E^(°^)^+°(*" /2 )- 

fc=0 

Now A; iterations of Proposition 14.31 shows that 



^(0, z) = [2(1 - *)] fc ^|(0, z) + £ /A* W^S=j(0, *), (4-5) 

i=i 

where the fj,k{z) are polynomials in z. With z = w + at 1 / 2 we have, by the product 
rule, 

y4(0, W + at^) = {n _^ k + j y (^ /2 r- 2k+ UM + o(t^). (4.6) 

We've also used the fact that q(w + at 1 / 2 ) = q(w) + o(l), because q has no zeros 
at to. From (|4.6[) we see that 

/ ilk («; + 0*1/2)^^(0, w + at 1 / 2 ) = o(r^), 
and can be ignored for j > 0. Putting this all together, we have 

£. f^ (o ' w+atl/2) - tn/2 1, w^W(- i ) fe [-(-™ i )] fc -"- 2 ^H+«(^ /2 )- 

Recalling the (physics) Hermite polynomial of order n, 

/ v v^, >t a™ _2fe n! 



fc!(n - 2fc)!' 
fc=0 v ; 



we obtain 



0(t, w + at 1 ' 2 ) = t n ' 2 [w(w - l)] n/2 H n (a/2^/w(w - l))q(w) + o(t n ' 2 ). 

H n is well known (e.g. 23, §3.3]) to have n distinct real roots, which, in the case 
that to < shows (|4.4j) . 

The case to = requires an additional argument. It will turn out, assuming 
4>(0,z) = z n q(z), that 4>(t,z) has n — 1 distinct negative zeros with distance of 
order t from zero, not t 1 ! 2 . 

To show this we again Taylor expand in t. By Proposition 14.31 we see that 



d k d d 2 



2 



d 

■z(l- z)—z n q(z 
az A 



where the operator z(l — z)d 2 is applied k times to z n q(z). If we evaluate this at 
z = at, we obtain 

^t(0,at) = [n(n - 1)] • • • [(n — fe + l)(n - fc)]M)"~V(0) + o(t n " fe ). 
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Thus the Taylor expansion gives 



«t,at)=f£J^a-» + o(f) 1 

where 



kl 

k=0 



[njfe = n(n — 1) • • • (n — k + 1) 

is the falling factorial. The sum can be rewritten as 

a n 7i!iFi[l-n,l,-l/a], 

where iJ*i [a, 6, x] is Kummer's hypergeometric function of the first kind. Now 
i-Fi[l — n, 1,— 1/x] has n—1 distinct negative zeros {zi,...z n -x} 22, pp. 103]. 
Hence for t small enough we can interlace n—1 zeros of <f>(t, z) between the points 
{0, Zi, . . . , z„_i}, and so the single remaining zero of (j>(t, z) near must also be 
real. □ 

Proof of "if" direction in Theorem \l.ll By Proposition ^. 13l the birth-death chain 
with constant birth and linear death rates preserves t-stability, and we just showed 
that the pure quadratic death chain preserves stability (and hence t-stability). 
However, the latter chain is no longer a Feller process, so we cannot immediately 
apply Trotter's product formula - as we did with reaction-diffusion processes and 
independent Markov chains - to combine the two processes. Indeed, it is well 
known that a pure quadratic death chain comes down from infinity in finite time, 
in the sense that liminffe^ooPt(fe, 1) > for each t > [TTj . 

We rectify this situation by considering the Banach space Z 1 (N) of absolutely 
summable sequences. Let X^, x[ 2 \ and be the birth-death chains with 

respective rates = &0j ^ — dik}, {(3^ — 0, 5^ = d 2 k(k — 1)}, and 

= b , 4 3) = d x k + d 2 k(k - 1)}. With 

p (i) (t)f(x) = J2 f(y) p ( x t l) - x\ x o ] = v) 

V 

as the (adjoint) strongly continuous contraction semigroups on Z 1 (N), we consider 
the infinitesimal generators as the I 1 limit 

nW/ = Um ^(*)/-/, 
tj.o t 

See P]5] for the theory of adjoint semigroups of Markov chains. 
Let 



D = {fel\N) 
D e = {fel\N) 
D(nW) ={/ e l\N) 



f(x) = for all but finitely many x}, 

\f(x)\ < Ce~ x }, C depending only on /, and 

\imr 1 (P {l) (t)f - /) exists as an I 1 limit.} 
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By explicit calculation, it can be seen that D C D e C £>(fiW) for each i, pW(i) : 
Do -> -D e , and for f E D e , 

= 40 i/(a . + 1} + 0W i/(a . _ i) _ + tfW]/^). 

By [9J Prop. 3.3 of Ch. 1], D e is a core for all three generators. Also, 

fiW + fiW = fiC3) on D e) 

so we can apply Trotter's product formula to conclude preservation of t-stability 
forX t (3) . □ 
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